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ABSTRACT 

Aims. We present a parameter study of the magnetohydrodynamical dynamo driven by cosmic rays in the interstellar medium (ISM) 
focusing on the efficiency of magnetic field amplification and the issue of energy equipartition between magnetic, kinetic and cosmic 
ray (CR) energies. 

Methods. We perform numerical CR-MHD simulations of the ISM using the extended version of ZEUS-3D code in the shearing box 
approximation and taking into account the presence of Ohmic resistivity, tidal forces and vertical disk gravity. CRs are supplied in ran- 
domly distributed supernova (SN) remnants and are described by the diffusion-advection equation, which incorporates an anisotropic 
diffusion tensor. 

Results. The azimuthal magnetic flux and total magnetic energy are amplified in majority of models depending on a particular choice 
of model parameters. We find that the most favorable conditions for magnetic field amplification correspond to magnetic diffusivity of 
the order of 3x 10 25 cm 2 s~', SN rates close to those observed in the Milky Way, periodic SN activity corresponding to spiral arms, and 
highly anisotropic and field-aligned CR diffusion. The rate of magnetic field amplification is relatively insensitive to the magnitude of 
SN rates in a rage of spanning 10% up to 100% of realistic values. The timescale of magnetic field amplification in the most favorable 
conditions is 150 Myr, at galactocentric radius equal to 5 kpc, which is close to the timescale of galactic rotation. The final magnetic 
field energies reached in the efficient amplification cases fluctuate near equipartition with the gas kinetic energy. In all models CR 
energy exceeds the equipartition values by a least an order of magnitude, in contrary to the commonly expected equipartition. We 
suggest that the excess of cosmic rays in numerical models can be attributed to the fact that the shearing-box does not permit cosmic 
rays to leave the system along the horizontal magnetic field, as it may be the case of real galaxies. 

Key words. Galaxies: ISM - magnetic fields - ISM: cosmic rays - magnetic fields - kinematics and dynamics - MHD 



1. Introduction 

An attractive idea of fast galactic dynamo has been proposed 
by Parker (1992). The idea relies on two ingredients: (1) cos- 
mic rays (CR) continuously supplied to the disk by supernova 
(SN) remnants and (2) fast magnetic reconnection which oper- 
ates in current sheets and allows to dissipate and relax the ran- 
dom magnetic field components in the limit of vanishing resis- 
tivity. Over the last decade we have investigated the different 
elements, physical properties and consequences of Parker's idea 
and scenario by means of analytical calculations and numeri- 
cal simulations (Hanasz & Lesch 1993, 1997, 1998, 2000, 2001, 
2003a, 2003b, Hanasz et al. 2002, 2004, 2006, Lesch & Hanasz 
2003; Otmianowska-Mazur et al. 2003, 2007; Kowal et al. 2003, 
2005) 

The first complete 3D numerical model of the CR-driven dy- 
namo has been demonstrated by Hanasz et al. (2004, 2006). In 
this paper we perform a parameter study of the CR-driven dy- 
namo model by examining the dependence of magnetic field am- 
plification on magnetic diffusivity, supernova rate determining 
the CR injection rate, temporal modulation of SN activity, grid 
resolution, and CR diffusion coefficients. 



The principle of action of the CR-driven dynamo is based on 
the cosmic ray energy supplied continuously by SN remnants. 
Due to the anisotropic diffusion of cosmic rays and the horizon- 
tal magnetic field configuration, cosmic rays tend to accumulate 
within the disc volume. However, the configuration stratified by 
vertical gravity is unstable with respect to the Parker instability. 
Buoyancy effects induce vertical and horizontal motions of the 
fluid and formation of undulated patterns - magnetic loops in 
frozen-in, predominantly horizontal magnetic fields. The pres- 
ence of rotation in galactic disks implies a coherent twisting 
of the loops by means of the Coriolis force, which leads to the 
generation of small-scale radial magnetic field components. The 
next phase is merging of small-scale loops by the magnetic re- 
connection process to form the large scale radial magnetic fields. 
Finally, the differential rotation stretches the radial magnetic 
field to amplify the large-scale azimuthal magnetic field com- 
ponent. The coupling of amplification processes of radial and 
azimuthal magnetic field components results in an exponential 
growth of the large scale magnetic field. The timescale of mag- 
netic field amplification, resulting from the action of the CR- 
driven dynamo, has been found (Hanasz et al. 2004, 2006) to be 
equal to 140- 250 Myr, depending on the galactocentric radius, 
which is comparable to the galactic rotation period. 
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The CR-dynamo experiments reported in the afore men- 
tioned papers rely on energy of CRs accelerated in SN remnants. 
Recently Gressel et al. (2008a,b) reported a series of non-ideal 
MHD simulations demonstrating dynamo action resulting from 
the SN-driven turbulence, in absence of CRs. Using a similar 
set of galactic disk parameters, with angular velocity enlarged 
4 times the value of 0.025Myr _1 , typical for the galactocen- 
tric radius of Sun, these authors found amplification of the large 
scale magnetic fields on a time scale of 250Myr . This indicates 
some similarity between the CR-driven dynamo and the dynamo 
driven by thermal energy output from supernovae. The similarity 
is presumably related to the buoyancy effect which can be com- 
monly attributed to the excess of both the thermal and cosmic 
ray energies in the disk volume. 

The magnitudes of galactic magnetic fields are usually es- 
timated from measurements of the radio synchrotron emission 
arising from the acceleration of cosmic ray electrons in the mag- 
netic field. To interpret the radio emission spectrum, it is usually 
assumed that the energy density in the magnetic field is of the 
same order of magnitude as the energy density in cosmic ray pro- 
tons (which are assumed to outnumber the electrons by 100:1, 
as they do in our Galaxy). There is however no compelling ev- 
idence of energy equipartition. Since the equipartition or min- 
imum energy assumption is one of the few ways to calculate 
radio source parameters, it is important to determine how rea- 
sonable the approach generally is. Recently Strong et al. (2007) 
and Snodin et al. (2005) raised again the question about the ap- 
plicability of the equipartition argument. 

From the observational point of view (Fitt and Alexander 
1993; Vallee 1995) the equipartition assumption seems to hold. 
In particular, Vallee's comparison of three different methods 
determining galactic magnetic field strengths (Faraday rotation 
method, equipartition method and Cosmic ray equipartition) 
shows that the equipartition fields are in a quite good agreement. 
On the other hand, Beck and Krause (2005) considered in detail 
a problem which was raised by Chi and Wolfendale (1993), that 
the commonly used classical equipartition or minimum-energy 
estimate of total magnetic field strengths from radio synchrotron 
intensities is of limited practical use because it is based on the 
hardly known ratio K of the total energies of cosmic ray protons 
and electrons and also has inherent problems. They present a re- 
vised formula using the number density ratio K for which they 
give estimates. For particle acceleration in strong shocks K is 
about 40 and increases with decreasing shock strength. Their re- 
vised estimate for the field strength gives larger values than the 
classical estimate for flat radio spectra with spectral indices of 
about 0.5-0.6, but smaller values for steep spectra and total fields 
stronger than about 10 //G. In very young supernova remnants, 
for example, the classical estimate may be too large by up to 10. 
On the other hand, if energy losses of cosmic ray electrons are 
important, K increases with particle energy and the equipartition 
field may be underestimated significantly. 

From a more global point of view the assumption of equipar- 
tition seems to be a very natural one. A thermodynamical system 
will always distribute the free energy to all degrees of freedom 
available, if the system has time to do so. In cases of acceler- 
ated particles diffusing in the large-scale and turbulent magnetic 
fields one can expect that at least the turbulent magnetic field 
(since it represents three degrees of freedom) is somehow virial- 
ized with respect to any other pressure term, like the cosmic ray 
pressure. This may not be true for the ordered magnetic fields, 
which are supposed to be amplified by the combined action of 
differentially rotating shear flows in the disk and some helical 
upward and downward motion driven either by cosmic rays pres- 



sure or any activity in the disk. In any case the connection of 
star formation activity accompanied by enhanced flux of cosmic 
rays and the amplification of large-scale magnetic fields, inher- 
ently raises the expectation that magnetic fields should not ex- 
hibit higher pressures than the cosmic rays. Moreover, one would 
rather expect to find magnetic fields whose pressure is somewhat 
lower than the pressure of the cosmic rays, if the cosmic rays 
present a source for the galactic dynamo. 

The paper is organized as follows: in Sect. [2] we describe the 
CR driven dynamo model and its numerical implementation, in 
Sect. [3] we present our simulation setup and describe parameters 
used in numerical simulations, in Sect. [4] we describe results, 
focusing on the effect of each parameter on magnetic field am- 
plification rate. We discuss the final saturated states of models in 
terms of equipartition between kinetic, magnetic and CR ener- 
gies. Finally, in Sect. [5] we conclude our paper. 

2. Description of the model 

Similarly as in the papers by Hanasz et al (2004, 2006) we take 
into account the following elements of the CR-driven dynamo: 

(1) The cosmic ray component, a relativistic gas, which is 
described by the diffusion-advection transport equation (see 
Hanasz & Lesch 2003b for the details of numerical algorithm). 
The typical values of the diffusion coefficient found from fitting 
to CR data (see eg. Strong et al 2007) are (3 -h5) x 10 28 cm V at 
energies ~ lGeV, and even larger values 10 29 cm 2 s~' are men- 
tioned (Jokipii 1999), however, we shall use reduced values in 
majority of simulations. 

(2) Following Giaccalone & Jokipii (1999) and Jokipii (1999) 
we presume that cosmic rays diffuse anisotropically along mag- 
netic field lines. The ratio of the perpendicular to parallel CR 
diffusion coefficients suggested by these authors is 5 %. 

(3) Localized sources of cosmic rays: supernova remnants ex- 
ploding randomly in the disk volume (see Hanasz & Lesch 
2004). We assume that each SN remnant supplies cosmic rays 
almost instantaneously, i.e. the comic ray input equal to 10 % of 
the canonical SN kinetic energy output (= 10 51 erg) for a single 
SN remnant is distributed over several subsequent time-steps. 

(4) Resistivity of the ISM (see Hanasz et al. 2002, Hanasz & 
Lesch 2003a, Tanuma et al. 2003) responsible for the onset of 
fast magnetic reconnection and topological evolution of mag- 
netic field lines. In this paper we apply the uniform resistivity 
and neglect the Ohmic heating of gas by resistive dissipation of 
magnetic fields. 

(5) Shearing boundary conditions and tidal forces following the 
prescription by Hawley, Gammie & Balbus (1995) aimed to 
model differentially rotating disks in the local approximation. 

(6) Realistic vertical disk gravity following the model of ISM in 
the Milky Way by Ferriere (1998). 

The set of equations describing the model of the CR-driven 
dynamo includes resistive MHD and the cosmic ray transport 
equations (see Hanasz et al., 2004): 



dp 
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— = V x (V x B) + i]AB, (4) 
at 

p = (y-l)e, r = 5/3 (5) 

where q = -d InQ/d InR is the shearing parameter, R is the dis- 
tance to galactic center, tj is the resistivity, y is the adiabatic 
index of thermal gas, the gradient of cosmic ray pressure Vp cr 
is included in the equation of motion (see eg. Berezinski et al. 
1990) and other symbols have their usual meaning. The uniform 
resistivity is included only in the induction equation (see Hanasz 
et al. 2002). The thermal gas component is currently treated as 
an adiabatic medium. 

The transport of the cosmic ray component is described by 
the diffusion-advection equation (see eg. Berezinski et al. 1990, 
Schlickeiser & Lerche 1985) 

+ V (e CT V) = V (RVe cr ) - Pcr (V ■ V) + Q SN , (6) 

where 2sn represents the source term for the cosmic ray energy 
density: the rate of production of cosmic rays injected locally in 
the SN remnants and 

Per = (Tcr - lKr, Tcr = 14/9. (7) 

The adiabatic index of the cosmic ray gas 7 cr and the formula for 
diffusion tensor 

K u = K ± 6 U + (K|| - KJmnj, m = Bi/B, (8) 
are adopted following the argumentation by Ryu et al. (2003). 

3. Numerical simulations 

3.1. Simulation setup 

In this paper we present a series of recent numerical simulations, 
whose aim is to search for the most favorable conditions for 
magnetic field amplification by means of the CR-driven dynamo. 
The presented numerical simulations were performed with the 
aid of Zeus-3D MHD code (Stone and Norman 1992 a,b) ex- 
tended with additions to the standard algorithm, that correspond 
to items (1) - (6) of Sect. [2] i.e. the cosmic ray component, 
treated as a fluid and described by the diffusion-advection equa- 
tion, including anisotropic CR diffusion tensor and cosmic ray 
sources - supernova remnants exploding randomly in the disk 
volume, resistivity of the ISM leading to magnetic reconnection, 
shearing-periodic boundary conditions, rotational pseudo-forces 
and a realistic vertical disk gravity. 

All simulations are performed in a Cartesian domain of size 
0.5 kpc x 1 kpc x 2 kpc in x, y, z coordinates, corresponding to 
radial, azimuthal and vertical directions, respectively. The basic 
resolution of the numerical grid is 50 x 100 x 400 grid cells in x, 
y and z directions, respectively, and for a smaller sample of simu- 
lations performed with larger values of CR diffusion coefficients 
the grid resolution is 25 x 50 x 200 grid cells. The boundary 
conditions are sheared-periodic in coordinate x, periodic in co- 
ordinate y and outflow on outer z-boundaries, with e cr — at the 
domain boundaries. The positions of SN are chosen randomly, 
with a uniform distribution in xy coordinates and Gaussian dis- 
tribution in z coordinate. 

The initial density distribution results from integration of the 
hydrostatic equilibrium equation, for the vertical gravity model 
of Ferriere (1998) and with the assumption of constant gas tem- 
perature across the disk, equal approximately 6000 K corre- 
sponding to the sound speed equal to about 7 km/s. The inte- 
gration procedure finds hydrostatic equilibrium for a given gas 
column density treated as an input parameter. 



The magnetic field strength, incorporated in the initial hy- 
drostatic equilibrium of gas, is defined through the parameter a 
denoting the ratio of initial magnetic to gas pressures. The initial 
cosmic ray pressure is equal to the initial gas pressure. 

The CR energy supplied to the system in SN remnants, ran- 
domly distributed around the disk midplane, implies that the CR 
pressure gradient force accelerates a vertical wind of thermal 
gas. To prevent significant mass losses from the computational 
domain, due to the vertical wind, we compensate the mass-loss 
Am after each timestep. The compensation mass is supplied as a 
mass source term, which is proportional to the initial mass dis- 
tribution 

Am 

Ap(x,y,z) = p (x,y,z), (9) 

w tot 

where m tot is the total gas mass in the computational domain 
and po(x, y , z) is the initial density distribution, resulting from the 
integration procedure of the hydrostatic equilibrium equation. 

3.2. Simulation parameters 

The basic input parameters, resulting from the assumed model 
are: the vertical gravity profile, local value of the galactic ro- 
tation and shear, gas column density and supernova rate. We 
adopt these parameters from the global model of ISM in Milky 
Way (Ferriere 1998) for the galactocentric radius R = 5kpc , 
where angular velocity is Q. — 0.05Myr _1 , gas column density 
2 = 27 x 10 20 cm ~ 2 and the realistic vertical gravity given by for- 
mula (36) in the aforementioned paper. The values of gas column 
density correspond in our simulations to the total density of all 
gas components in Ferriere (1998), while SN-rate is the rate of 
type II supernovae. We assume for simplicity that all SN explo- 
sions appear as single supernovae, and that vertical distribution 
of SN explosions is Gaussian, with a fixed half-width equal to 
100 pc. 

In addition to the mentioned relatively well established lo- 
cal disk parameters, there is a group of less known quantities 
like effective magnetic diffusivity, CR diffusion coefficients and 
efficiency of conversion of SN kinetic energy into cosmic ray 
energy. We assume the standard 10% value of kinetic CR en- 
ergy conversion efficiency and the magnetic diffusivity and CR 
diffusion coefficients varying in a wide range. 

In this paper we present the results of five simulation series 
A-E. The summary of all variable simulations parameters for the 
whole set of simulations is presented in Table Q] 

In simulation series A (runs Al - A5) we examine effects 
of magnetic diffusivity variations on magnetic field amplifica- 
tion by applying rj in the range H- 10 3 pc 2 Myr _I corresponding 
to t 3 x 10 26 cm 2 s~' in CGS units. We define the magnetic 
Reynolds number Rm = LlQ/j] for reference, as in Gressel et 
al. (2008b), where L y = lOOOpc is the domain size in y direc- 
tion. Moreover we assume continuous and time-invariable sup- 
ply of CRs in SN remnants. The simulation runs Al, A2, A3 
are the same as the runs B, C, and D, respectively, discussed by 
Otmianowska-Mazur et al (2007). We note, for comparison, that 
the commonly adopted value of turbulent diffusivity in the ISM 
is 77turb - l/3v turb L turb ~ 10 26 cm V 1 1/3 X 10 3 pc 2 Myr -') 
for v tur b = 10km s _1 and L tur b = lOOpc . We note that the adopted 
values of magnetic diffusivity exceed the value corresponding to 
the Spitzer resistivity (^ 10 8 cm 2 s~', see Parker 1992) by 15- 
18 orders of magnitude. The relative smallness of the Spitzer 
resistivity implies that an anomalous resistivity, considered as 
a subscale phenomenon, has to be invoked in order to explain 
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Table 1. Parameters of simulations presented in this paper. Subsequent columns show: simulation name, initial ratio of magnetic to 
gas pressure a, magnetic diffusivity 77, magnetic Reynolds number Rm, surface frequency of SN explosions, presence of SN-rate 
modulation, parallel K\\ and perpendicular K ± CR diffusion coefficients, and grid resolution in x, y and z directions. 



dissipation of the small scale magnetic fluctuations in the ISM. 
Following Parker (1992), we assume that reconnection rates in 
the ISM are comparable to those predicted by the Petscheck's 
fast reconnection model, i.e. the magnetic cutting speeds are of 
the order of va / log(Rivi) rather than va I VRvi typical for the slow 
Parker-Sweet reconnection model, where Rm is the Lundquist 
number or the magnetic Reynolds number. 

In the simulation series B (runs Bl - B5) we apply the same 
range of magnetic diffusivity values, but CRs supply is modu- 
lated in a manner mimicking passages of subsequent spiral arms, 
regulating the star formation rate and subsequently SN-rate. The 
effect of spiral arms is modelled (see Hanasz et al. 2006) by sup- 
plying cosmic rays in intermittent periods of 25 Myr, with SN 
rate equal to 4x the reference /sn, and followed by periods of 75 
Myr without any SN activity. The time-averaged supernova rate 
is in this case equal to the reference /sn- 

In the simulation series C (runs C1-C5) we apply a constant 
magnetic diffusivity rj = lOOpc 2 Myr _1 and vary the surface fre- 
quency of SN explosions in the range of 15 - 500 kpc ~ 2 Myr _1 , 
assuming modulated CR supply as in the simulation series B. In 
the simulation series D (runs Dl - D2) we repeat simulations A4 
and B4, respectively, with the grid resolution reduced twice in 
each spatial direction. 

Due to the CFL timestep limitation of the currently used 
explicit algorithm of CR diffusion, the applied values of CR 
diffusion coefficients are scaled down with respect to the re- 
alistic values. The timestep limitation ensuring stability of ex- 
plicit numerical schemes applied to the diffusion equation is 
At < 0.5(Ax) 2 /K, where K is the diffusion coefficient. The 
timestep becomes prohibitively short when the diffusion coef- 
ficient is very large or the spatial step is too small. For this 
reason the CR diffusion coefficient has been reduced in simu- 
lation series A-D, by about one order of magnitudes, with re- 
spect to the mentioned realistic values 3 r 6 X 10 28 cm 2 s~' ^ 
1 -r 2 ■ 10 5 pc 2 Myr _1 . The fiducial values of the parallel and 
perpendicular diffusion coefficients applied in simulation series 



A-D are respectively: K\\ = 1 X 10 4 pc 2 Myr 1 =s 3 x 10 27 cm 2 s 1 
and K ± = 1 x 10 3 - 3 x 10 26 cm V 1 . 

Finally, in the simulation series E (runs El - E6) we increase 
the parallel K\\ and perpendicular K± CR diffusion coefficients 
by factors 3 and 10, with respect to the fiducial values, in or- 
der to examine magnetic field amplification for more realistic 
magnitudes of these quantities. In this way we apply realistic 
CR diffusion coefficients in a few single simulation runs, yet the 
maximum CR diffusion coefficients used do not reach the upper 
range of realistic values, of the order of 10 29 cm 2 s~', mentioned 
in the literature. 

4. Results 

In Fig. [TJ we show the distribution of cosmic ray gas together 
with magnetic field vectors (panel (a)), and thermal gas density 
together with gas velocity vectors (panel (b)) in the yz-slice taken 
forx = 0atf= 1000 Myr. 

One can notice in panel (a) that the dominating horizontal 
magnetic field component in the disk volume is undulated in a 
manner resembling the effects of Parker instability. The cosmic 
ray energy density is well smoothed by the diffusive transport 
in the computational volume. The vertical gradient of the cos- 
mic ray energy density is maintained by the supply of cosmic 
rays around the equatorial plane in the disk in the presence of 
vertical gravity. The cosmic ray energy density is expressed in 
units in which the thermal gas energy density corresponding to 
p = 1cm ~ 3 and the isothermal sound speed c s = 7km s _1 is 
equal to 1 . The velocity field together with gas density is shown 
in panel (b). It is apparent that the distribution of gas is signifi- 
cantly less smooth that the distribution of cosmic rays. 

In order to examine the structure of the large-scale field we 
show in panel (c) the horizontally averaged magnetic field com- 
ponents (B x (z)) and and (B y (z)). A striking property of the mean 
magnetic field configuration is the almost exact coincidence of 
peaks of the oppositely directed radial and azimuthal field com- 
ponents. This feature resembles to the standard picture of an aa>- 
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dynamo: the azimuthal mean magnetic component is generated 
from the radial one and vice versa (see Lesch and Hanasz (2003) 
for a corresponding simple analytical model). 

In panel (d) we show the horizontally averaged vertical ve- 
locity component (V z ) and its fluctuations (SV Z ). It is apparent 
that the bulk speeds of the wind, driven by the vertical gradient of 
CR pressure, reach 65 ± 15km /s at z = ±2kpc . Vertical system- 
atic winds blowing with bulk speeds, comparable to rotational 
galactic velocities, influence large-scale structures of galactic 
magnetic fields and are observed in external starburst galaxies 
like NGC 253 (see Beck 2007, Hessen et al. 2008). 

By varying the parameters discussed in the previous section 
we intend to determine the regions of parameter space in which 
the magnetic field amplification is the most efficient. The am- 
plification of the regular magnetic field is identified with the 
amplification of the total magnetic energy in the computational 
domain, associated with the amplification of the azimuthal mag- 
netic flux. The magnetic flux shown in subsequent plots repre- 
sents an azimuthal flux averaged over all xz slices (spread in y 
direction) through the discretized computational domain. In the 
subsequent subsections we present the parameter study of the 
CR-driven dynamo, focusing on the efficiency of magnetic field 
amplification and the issue of equipartition between magnetic, 
kinetic, and CR energies. 

4.1. Dependence of magnetic field amplification on magnetic 
diffusivity 

As a first step in our parameter study of the cosmic-ray driven 
dynamo we examine, in simulation series A, the effect of mag- 
netic diffusivity on the efficiency of magnetic field amplification. 
Time evolution of magnetic energy and magnetic flux are shown 
in Fig. [2] Magnetic flux plotted in the left panel of Fig.[2]is scaled 
in the following way. The initial magnetic field induction is de- 
fined by the parameter a = p mag l p gas , shown in the second col- 
umn of Table. [1] where we apply Parker's convention to assign 
the inverse of plasma beta as a. The adopted values of a are 10~ 4 
and 10~ 2 in different simulations, while a - I means magnetic 
pressure equal to the thermal gas pressure. We scale magnetic 
flux in such a way that a = 1 corresponds to the azimuthal mag- 
netic flux (D fl = 1 . The total magnetic energy plotted in the right 
panel of Fig. [2]is scaled with respect to the time-averaged total 
kinetic energy in the computational domain. The latter quantity 
appears to fluctuate around a mean value, which is practically 
time-invariant for all simulation runs, thus we find this kind of 
scaling convenient. The scaling described above will be applied 
to all subsequent plots of magnetic flux and magnetic energy. 

The evolution of magnetic energy and magnetic flux in mod- 
els of simulation series A, represented by different curves in 
Fig. |2] demonstrate that magnetic field amplification strongly 
depends on magnetic diffusivity. In the case of vanishing ex- 
plicit resistivity, magnetic energy grows up by about 2.5 orders 
of magnitude during the first 500 Myr, but magnetic flux is am- 
plified only by a factor of 3 during this period and it fades later 
on. This effect can be interpreted as a predominant growth of 
the small-scale turbulent magnetic field component with a lit- 
tle contribution of large-scale magnetic field amplification (see 
Otmianowska-Mazur et al. 2007 for a more extended analysis of 
the simulations presented in this section). One should remember, 
however, that numerical resistivity, always present in numerical 
MHD simulations, may influence to some extent the behavior 
of the simulation run A, corresponding to r\ = 0. The amount 
of numerical magnetic diffusivity 77 ^ 0.7pc 2 Myr ~' has been 



quantified, for the present grid resolution, on the base of Parker 
instability simulations by Kowal et al (2003). 

When magnetic diffusivity is increased up to 77 = 
100 pc 2 Myr _1 the efficiency of magnetic field amplification 
increases. For 77 = lOOpc 2 Myr ~' (run A4) the growth of mag- 
netic flux persists until t = 1000 Myr and saturates thereafter. 
For smaller values of 77 the growth rate is smaller and the maxi- 
mum values of magnetic flux are smaller than those attained for 
77 = lOOpc 2 Myr _1 . In cases of vanishing or small explicit dif- 
fusivity (runs Al and A2), magnetic energy grows up initially 
faster than in the case of larger resistivity (runs A3 and A4). 
This behavior means that low resistivity enables initially faster 
growth of the random magnetic field component, while for larger 
resistivity random magnetic fields are quickly dissipated. The 
growth of the total magnetic energy follows closely, in the lat- 
ter case, the growth of the mean magnetic flux. It is also appar- 
ent that amplification of magnetic flux and magnetic energy for 
77 = lOOOpc 2 Myr _1 (run A5) is significantly reduced with re- 
spect to the other runs. 

In order to explain the physical mechanism that controls 
magnetic field amplification through the magnitude of magnetic 
diffusivity we plot the ratio of total (volume integrated) energies 
of vertical to azimuthal magnetic field components in Fig. [3] 
It is apparent that the energy of vertical magnetic field dom- 
inates in all lower magnetic diffusivity runs A1-A3 (77 = 0, 1 
and lOpc 2 Myr _1 ) and is comparable to the energy of azimuthal 
magnetic field in run A4 (77 = lOOpc 2 Myr _1 ) providing the 
strongest magnetic field amplification. Among all simulations of 
series A only simulation A4 reaches energetic equipartition be- 
tween magnetic field and gas kinetic energy as a result of the 
amplification process. 

On the other hand, in the case of large magnetic diffusivity 
(run A5, 77 = lOOOpc 2 Myr _1 ) energy of the vertical magnetic 
field component remains much smaller than energy of the az- 
imuthal magnetic field and magnetic field amplification does not 
occur. This fact can be interpreted by resistive damping of the 
undulatory mode of the Parker instability in favor of the inter- 
change mode, which does not contribute to the dynamo action. 
The above finding indicates that the most favorable conditions 
for magnetic field amplification corresponds to approximately 
equal energies in vertical and azimuthal magnetic fields in the 
case of buoyancy driven dynamo. 

To demonstrate the effect of resistivity on kinetic and mag- 
netic turbulent spectra we compute the Fourier transforms of 
kinetic and magnetic energy densities, as in the paper by 
Otmianowska-Mazur et al (2007). The results are shown in Fig|4] 
The highly anisotropic nature of turbulence is reflected by differ- 
ent lines representing Fourier transforms in x, y and z-directions. 
We find that the kinetic spectra, which are generally close to 
the Kolmogorov spectrum oc k~ 5/3 , depend rather weakly on 
magnetic diffusivity. The large values of magnetic diffusivity 
77 = 100 t lOOOpc 2 Myr _I lead only to the damping of the 
short-wavelength components of the Fourier spectrum computed 
in y-direction. 

The magnetic spectra appear to be much more sensitive to 
variations in magnetic diffusivity. The plots obtained for runs Al 
and A2 exhibit practically identical spectra in all directions. This 
means that the diffusivity of lpc 2 Myr _1 does not change the 
results with vanishing explicit resistivity, or in other words the 
numerical resistivity of the code corresponds to the resistivity of 
Run A2. The effect of resistivity is apparent through the reduced 
amplitude of large-/c modes for 77 = lOpc 2 Myr _1 (Run A3) and 
an apparent cutoff in magnetic spectra around k/2n 0.2pc~ l 
for 77 = lOOpc 2 Myr _1 . Further increase of magnetic diffusivity 
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Fig. 1. Exemplary plots illustrating the state of the system at t — lOOOMyr for simulation A4. In the first two panels we present slices 
through the computational volume in the yz-plane for x = 0. Panel (a) shows cosmic ray energy density with vectors of magnetic 
field, panel (b) shows gas density with velocity vectors. In panel (c) we plot horizontally averaged x and y components of magnetic 
field, and in panel (d) horizontally averaged vertical velocity component and its fluctuations. 



up to j] — lOOOpc 2 Myr ~' leads to the signifficant reduction of 
amplitudes of all modes for the Fourier tranforms performed in 
z-direction, steppening of the whole spectrum in z-direction and 
a surprising effect of flattening of the spectrum in x-direction. 
The latter effect may indicate a qualitative change of the phys- 
ical nature of the modes, which can plausibly attributed to the 
mentionned enhancement of the exchange mode of Parker insta- 
bility. 

We briefly note that we neglect any small scale dynamics 
of the helical MHD turbulence. This assumption allows the uni- 
formity of our diffusion coefficients, like the magnetic diffusiv- 
ity. Detailed investigations of the influence of small scale helical 
MHD turbulenceon galactic dynamos are presented by Maron & 
Blackman (2002) and Maron, Cowley and McWilliams (2004). 

One can further interpret these results in terms of topolog- 
ical evolution of magnetic field which is controlled by resistiv- 
ity. The topology of magnetic field lines determines the paths of 
anisotropic cosmic ray transport. For low values of resistivity the 
buoyancy of cosmic rays leads to opening of magnetic field lines 
through upper- and lower-z boundaries. This implies that diffu- 
sive escape of cosmic rays, along the open magnetic field lines, 
dominates over the buoyancy and limits the effect of Coriolis 
force, which is responsible for dynamo action. 

4.2. The effect of spiral arms 

In this section we describe the simulation series B performed for 
the same set of simulation parameters as for series A with an ex- 



ception that currently the cosmic ray energy input is modulated 
in time by a step function. Motivation for this kind of CR supply 
is the presence of spiral arms in disk galaxies. We assume that 
SNe explode in arms with the rate, which is proportional to star 
formation rate. We assume that arms pass trough the volume of 
our local computational domain once per lOOMyr and that the 
arm passage takes 25Myr , i.e. starting at t — we supply CRs 
for the first 25Myr of the lOOMyr and stop CR supply for the 
remaining 75Myr . We enhance the SN rate in arms 4 times, so 
that the average SN rate over the whole period of density wave 
remains the same as in the simulation series A. The evolution of 
mean magnetic flux and energy is presented in Fig. [5] 

We find that in the present set of simulations, both magnetic 
flux and magnetic energy grow up faster than in the case of sim- 
ulations without CR modulation. Magnetic field amplification is 
now apparent even in simulations with r] = 0. The only exception 
is the j] = lOOOpc 2 Myr _1 simulation which does not show a no- 
ticeable amplification of magnetic field, even in the presence of 
CR modulation. This indicates that the temporal modulation of 
SN rate acts in the same way as increasing magnetic diffusivity 
in the range of rj = 4- lOOpc 2 Myr _I . 

To interpret the above results we suggest the following sce- 
nario: CRs supplied to the system trigger Parker instability and 
leave the disk volume via combined buoyant and diffusive trans- 
port. Vertical magnetic loops form efficiently during the period 
of enhanced SN activity, but later on, in absence of CR pertur- 
bations, magnetic field tends to relax before the next spiral arm 
passage. In absence of CR forcing, in the inter-arm regions even 
a small resistivity is sufficient to relax magnetic field structure 
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Fig. 2. Time evolution of azimuthal magnetic flux and total magnetic energy for different values of magnetic diffusivity in simulation 
series A. The curves represent respectively cases of 77 = (Al), 77 = 1 (A2), 77 = 10 (A3), 77 = 100 (A4) and 77 = 1000 (A5) in units 
pc 2 Myr -1 . 




Fig. 5. Time evolution of azimuthal magnetic flux and total magnetic energy for different values of magnetic diffusivity, in pres- 
ence of temporal modulations of SN-rate mimicking the presence of spiral arms in the simulation series B. The curves represent 
respectively cases of 77 = (Bl), 77 = 1 (B2), 77 = 10 (B3), 77 = 100 (B4) and 77 = 1000 (B5) in units pc 2 Myr _1 . 



to the horizontal more regular configuration which suppresses 
excessive losses of CRs via the diffusive transport. Thus the ef- 
ficiency of magnetic field amplification is enhanced. 

4.3. Dependence of magnetic field amplification on SN-rate 

In all runs of the simulation series A and B we adopted 
the fiducial value of /sn - 130kpc ~ 2 Myr _I derived from 
the global model of Milky Way by Ferriere (1998) at galac- 
tocentric radius Rc = 5kpc . In this section we describe 
simulation series C performed for different SN rates /sn = 



15,30,60,250 and 500kpc ~ 2 Myr _1 , together with run B4 
(fsN - 130kpc ~ 2 Myr _1 ), to examine the effect of SN rate on 
magnetic field amplification. The results are shown in Fig. [6] In 
all cases the SN input is modulated in a manner described in 
Sect. El 

One can notice that the magnetic field amplification rates and 
the final saturation levels of both magnetic flux and magnetic en- 
ergy grow with fs n as long as the SN rate is lower or equal to the 
fiducial, realistic value of fsN = 130kpc ~ 2 Myr _1 at galactocen- 
tric radius Rq = 5kpc , when all other disk parameters are fixed. 
The e-folding times of magnetic flux deduced from the left panel 
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Fig. 6. Time evolution of azimuthal magnetic flux and total magnetic energy for different values of SN rate applied in simulation 
series C, together with run B4, in the presence of temporal modulations of SN-rate. Line assignments are respectively: fsN = 15 
(CI), f SN = 30 (C2), f SN = 60 (C3), f SN = 130 (B4), f SN = 250 (C4), and f SN = 500 (C5) supernova explosions per squared kpc 
per Myr. 
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those simulation runs of series C for which the amplification 
holds. 

When the SN rate is doubled, with respect to the fiducial 
value, then only a short period of magnetic field amplification 
is observed, until t — 500 Myr , and if SN rate is doubled once 
again then magnetic field decays. The above results show that 
magnetic field amplification holds in a wide range of SN rates 
and that the realistic values of SN rates are optimal for the galac- 
tic dynamo process. Similarly as in Sect. I4.1l we show the ratio 
of energies in vertical to horizontal magnetic field components in 
Fig. H We find that for SN rates up to f SN = 130kpc ~ 2 Myr 1 
the efficient magnetic field amplification is associated with the 
ratio of vertical to horizontal magnetic field energies fluctuat- 
ing around one, and in the case of excessive CR supply (f$ n = 
250kpc ~ 2 Myr ~' and more) the vertical magnetic field energy 
dominates and magnetic field ceases to grow. 



Fig. 3. Time evolution of the ratio of energies of vertical to hor- 
izontal magnetic field components for different values of mag- 
netic diffusivity in the simulation series A. Line assignments are 
the same as in Fig. [2] 



of Fig. |6]are respectively 150 Myr for f SN = 130kpc 2 Myr 1 
(Run B4) and 190 Myr for f SN = 15kpc ~ 2 Myr _1 (Run CI). We 
note therefore that the magnetic field amplification rate is rela- 
tively insensitive to the magnitude of SN rate within the range 
of SN rates spanning one decade below the fiducial value. We 
note also that magnetic field amplification saturates at the level 
of equipartition of magnetic and kinetic energies in the case of 



4.4. Dependence of magnetic field amplification the grid 
resolution 



In order to check the influence of the grid resolution on simula- 
tion results we increase the cell size to (20pc ) 3 in simulations 
Dl and D2 and apply the same parameters as in simulations A4 
and B4, respectively. In Fig. [8] we show the evolution of the to- 
tal flux of the azimuthal magnetic field component and the total 
magnetic energy for simulations Dl and D2 together with anal- 
ogous curves for simulations A4 and B4, shown previously in 
Figs. |2]and|5] It is apparent that the results obtained at both res- 
olutions are very similar, although a slightly faster growth of 
magnetic field is observed in simulations performed at the lower 
resolution, which can be explained by somewhat larger numeri- 
cal resistivity. 
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Fig. 8. Time evolution of the azimuthal magnetic flux and the total magnetic energy for simulations with grid resolutions (lOpc ) 3 
(runs A4 and B4) and (20pc ) 3 (runs Dl and D2). 



4.5. Dependence of magnetic field amplification on CR 
diffusion coefficients 

The aim of simulation series E is to examine the effect of vari- 
ations of the CR diffusion coefficients on magnetic field am- 
plification. All simulations of series E are performed with the 
resolution (20pc ) 3 . The reduced grid resolution makes it possi- 
ble to enlarge CR diffusion coefficients to realistic values, while 
preserving acceptable timesteps for the explicit integration algo- 
rithm of the CR diffussion-advection equation. 

In the present simulation series E we vary parallel and per- 
pendicular diffusion coefficients choosing different pairs from 
the set: K ]} = 1 x 10 4 , 3 x 10 4 and 1 x 10 5 pc 2 Myr "\ and 
K ± = lx 10 3 , 3 x 10 3 and 1 x 10 4 pc 2 Myr _1 . The results of new 
simulations as compared to the run Dl are presented in Fig. [9] 

The results of simulation series E can be summarized as 
follows. We note, that magnetic field growth rate and the sat- 
uration values of magnetic flux and energy depend particularly 
on the choice of K\\ and K ± . When K\\ is increased 3 and 10 
times with respect to Dl, the initial growth of magnetic field 
becomes slightly faster and the saturation level becomes lower 
by a factor of 2 4- 3, provided that K ± is not too large. For K\\ = 
3xl0 4 pc 2 Myr~' amplification holds for K ± = lxl0 3 pc 2 Myr _I 
and 3 x 10 3 pc 2 Myr _1 , but for K ± = 1 x 10 4 magnetic field de- 
cays. Similarly, for K\\ = 1 X 10 5 pc 2 Myr _1 amplification holds 
for K ± = 1 x 10 3 pc 2 Myr _1 and K ± = 3 x 10 3 pc 2 Myr _1 , 
but for K ± = 1 X 10 4 pc 2 Myr _1 we find only initial growth 
until t = 450 Myr and decay thereafter. The present results 
indicate that magnetic field amplification is possible only for 
K ± < 3xl0 3 pc 2 Myr~' 10 2 ^ 'cm 2 s' 1 .Therefore the anisotropy 
ofCR diffusion seems to be a crucial condition for magnetic field 
amplification in the process of the CR-driven dynamo. 

In the subsequent Fig. [10] we show the energy ratio of ver- 
tical to azimuthal magnetic field components for runs E4, E5 
and E6, corresponding to three different values of K ± and K\\ = 
1 x 10 5 pc 2 Myr _1 . Comparing Figs.|9]and[10]we find that in the 
case of two simulation runs E4 and E5 (two smaller values of 
K ± ) the energy ratio of vertical to azimuthal magnetic field com- 



ponents varies in the range ~ 0.3 -r 2, corresponding to efficient 
growth magnetic energy. For the largest value of K ± = 10 4 the 
magnetic energy ratio becomes occasionally larger by an order 
of magnitude, magnetic field decays. 

4.6. The issue of energy equipartition 

The results presented so far demonstrate that magnetic fields 
amplified by the CR-driven dynamo saturate near equipartition 
of magnetic and gas kinetic energies. It is commonly expected 
that CRs remain in energetic equipartition with gas and magnetic 
field as well. 

To investigate the relation of CR to other forms of energies, 
we plot in Fig. [TT]the time evolution of the ratio of CR to the 
time-averaged turbulent kinetic energy (we subtract the kinetic 
energy of the large-scale shear flow) for different values of CR 
diffusion coefficients (the simulation series E). We find that, de- 
pending on diffusion coefficients, CR energy is larger than the 
turbulent kinetic energy by a factor of 10 - 50, while the mag- 
netic field energy, according to the results presented in Fig. [9] is 
close to the gas turbulent energy. 

It is apparent that for the first few hundred Myr CRs accu- 
mulate quickly in the disk since they are trapped by horizontal 
magnetic field. The ratio of CR to kinetic energies saturates as 
soon as vertical magnetic field component appears, due to buoy- 
ancy, enabling diffusive transport of CRs out of the disk. 

When the cosmic ray diffusion coefficients are larger, one 
can find that the ratio of CR to kinetic energies is lower. We note, 
that cosmic rays find an easier way to leave the disk, when the 
parallel and perpendicular diffusion coefficients are larger. The 
results shown in Fig. QT]indicate that the ratio of CR to kinetic 
energies anticorrelates with both: the parallel and perpendicular 
CR diffusion coefficients. Due to the mentioned timestep limi- 
tation, in the simulations presented in this paper we could only 
adopt the values of K« reaching at most 3 x 10 28 cm 2 s _1 , that are 
still smaller than 10 cm 2 s~' mentioned by other authors (e.g. 
Jokipii 1999). Although the tendency of lowering the CR en- 
ergy with the magnitude of the parallel CR diffusion coefficient 
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Fig. 9. Time evolution of the azimuthal magnetic flux and the total magnetic energy for different values of the parallel and perpen- 
dicular CR diffusion coefficients. Thin lines are used for K\\ = 1 x 10 4 (runs Dl and El), mid lines are used for K\\ — 3 x 10 4 (runs 
E2 E3 and E4) and thick lines are used for K\\ = 1 x 10 s (runs E5 E6 and E7). Full lines denote K ± = 1 x 10 3 (runs Dl, E2 and E5), 
dotted lines denote K ± = 3 x 10 3 (runs El, E3 and E6), dashed lines Kj_ = 1 x 10 4 (runs E4 and E7). All diffusion coefficients given 
in units pc 2 Myr _1 . 



seems promising, one should not expect that values as large as 
Aj| =s 10 29 cm 2 s _1 will reduce the problem of CR energy excess. 

Another factor, which may significantly influence the rela- 
tion between CR and other form of energies, is the choice of pe- 
riodic and shear-periodic boundary conditions in the horizontal 
directions of the computational box. In real galactic disks, CR 
diffusion is expected along horizontal magnetic field lines. In 
the case of periodic-type boundary conditions, CRs are trapped 
in the disk volume by predominantly horizontal magnetic field. 
This kind of trapping can be released, however, only in the global 
galactic disk simulations. 



5. Summary and conclusions 

In the present paper we described an extensive series of simu- 
lations and presented parameter study of the CR-driven dynamo 
in a galaxy, characterized by the essential parameters typical for 
the Milky Way at galactocentric radius R c = 5kpc . In the pre- 
sented study the magnetic diffusivity, as well as the parallel and 
perpendicular CR diffusion coefficients have been considered as 
free parameters, and dedicated simulation series have been per- 
formed to investigate their influence on the efficiency of the CR- 
driven dynamo process. The results of the parameter study can 
be summarized as follows: 

(1) The magnitude of magnetic diffusivity influences essen- 
tially the efficiency of magnetic field amplification. The most 
favorable value of magnetic diffusivity is lOOpc 2 Myr _1 

3 x 10 25 cm 2 s~', the value comparable, although lower than the 
the value of turbulent diffusivity of ISM deduced from observa- 
tional data. 

(2) The efficiency of magnetic field amplification is en- 
hanced by temporal modulation of the CR supply. An effect of 
this kind may be associated with the periodicity of star formation 
and supernova activity induced by spiral arms. The enhancement 



is apparent at lower values of magnetic diffusivity and is less sig- 
nificant at the optimal value of r/ = lOOpc 2 Myr _1 . 

(3) Magnetic field amplification rate is relatively insensitive 
to the magnitude of SN rate within the range of SN rates span- 
ning one decade below the value /sn = 130kpc ~ 2 Myr _1 typical 
for galactocentric radius R G = 5kpc . We note also that mag- 
netic field amplification saturates at the level of equipartition of 
magnetic and kinetic energies for all supernova rates for which 
amplification holds. Magnetic field is no longer amplified, if SN 
rate is further enhanced by factors 2 and 4 with respect to the re- 
alistic value, while other quantities (like e.g. gas column density) 
remain fixed. 

(4) Magnetic field amplification in the CR-driven dynamo 
relies on anisotropic diffusion of cosmic rays. From the limited 
set of simulations of series E, one can deduce that the magnetic 
field amplification is possible only for K ± < 3 x 10 3 pc 2 Myr _1 
10 27 cm 2 s~ l , and for all considered values of the parallel diffu- 
sion coefficient K\\ in the range 3 4- 30 x 10 27 cm 2 s _1 . Therefore, 
the 5% ratio of the perpendicular to parallel diffusion coefficients 
postulated by Giaccalone & Jokipii (1998) falls within the am- 
plification range. 

(5) By varying various parameters: magnetic diffusivity, su- 
pernova rate and the CR diffusion coefficients, we have found 
that the favorable conditions for magnetic field amplification 
correspond to approximately equal energies of the vertical and 
azimuthal magnetic field components in the case of buoyancy 
driven dynamo. An excess or deficit of vertical magnetic field 
with respect to the azimuthal one corresponds to significantly 
less efficient amplification or even decay of magnetic field. 

(6) We note the problem indicated previously by Snodin et 
al (2005), that in all simulations the CR energy in the computa- 
tional domain exceeds the turbulent kinetic energy and magnetic 
energy by more than one order of magnitudes. Moreover, the 
lowest ratios of CR to kinetic energies emerge for the largest 
values of the parallel diffusion coefficient. It seems not plausi- 
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Fig. 7. Time evolution of the ratio of energies of vertical to hor- 
izontal magnetic field components for different (modulated) SN 
rates f SN = 15 (CI), f SN = 130 (B4) and f SN = 250 (C4) super- 
nova explosions per squared kpc per Myr. 
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Fig. 4. Kinetic (left column) and magnetic (right column) spectra 
computed for t = lOOOMyr (runs A1-A4) and t = 700Myr (Run 
A5), separetelly in x, y and z-directions (full, dotted and dashed 
thin lines, respectively). Lines representing the k~ 5/3 slope (thick 
full lines) are shown for comparison. 



Fig. 10. Time evolution of the ratio of energies of vertical to 
horizontal magnetic field components for different values of 
the perpendicular CR diffusion coefficients K\\ = 3 x 10 4 and 
^ = 1X 10 4 (E4), K\\ = 1 X 10 5 and K ± = 1 X 10 3 (E5) and 
K\\ = 1 X 10 5 , K ± = 3 X 10 3 (E6). All diffusion coefficients Gaven 
in units pc 2 Myr _1 . 



ble, however, that an enlargement of diffusion coefficients up to 
fully realistic values will reduce the excess of cosmic rays en- 
ergy in the disk. It seems also, that the ratio of CR to other forms 
of energy in the ISM is not yet well restricted on observational 
grounds (Strong et al 2007). On the other hand, the currently 
used shearing box approximation does not permit CRs to leave 
the disk by means of diffusion along the predominantly horizon- 
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Fig. 11. Time evolution of the ratio of CR to time-averaged ki- 
netic energy for different values CR diffusion coefficients in sim- 
ulation series E. Line assignments are the same as in Fig. [9] 



tal magnetic field. Therefore, we suggest that subsequent work 
on the CR-driven dynamo, aimed to solve this problem, should 
be done in the framework of global galactic disk simulations. 

Acknowledgements. This work was supported from the Polish Committee for 
Scientific Research (KBN) through the grants PB 0656/P03D/2004/26 and 
2693/H03/2006/31. 



References 

Beck R., 2007. larXiv:071 1.47001 

Beck, R, Krause, M., 2005, Astronomische Nachrichten 326 

Berezinskii, VS., Bulanov, S.V., Dogiel, V.A., Ginzburg, V.L., Ptuskin, VS., 

Astrophysics of cosmic rays, Amsterdam: North-Holland, 1990. 
Chi, X., Wolfendale, A.W., Nature 362, 610 (1993) 
Ferriere, K. 1998, ApJ, 497, 759 
Fitt, A.J., Alexander, P., MNRAS, 261, 445 (1993) 
Giacalone, J., Jokipii, R.J., 1999, ApJ, 520, 204 
Gressel, O., Ziegler, U., Elstner, D., R "udiger, G. 2008, AN, 329, 61 
Gressel, O., Elstner, D., Ziegler, U., Rtidiger, G. 2008, A&A, 486L, 35 
Hanasz, M., Lesch, H. 1993, A&A, 278, 561 
Hanasz, M., Lesch, H. 1997, A&A, 321, 1007 
Hanasz, M., Lesch, H. 1998, A&A, 332, 77 
Hanasz, M. Lesch, H. 2000, ApJ, 543, 235 
Hanasz, M. Lesch, H. 2001, Space Sci. Rev., 99, 231. 
Hanasz, M., Otmianowska-Mazur, K., Lesch, H. 2002, A&A, 386, 347 
Hanasz, M., Lesch, H. 2003a, A&A, 404, 389 
Hanasz, M., Lesch, H. 2003b, A&A, 412, 331 

Hanasz, M., Kowal, G., Otmianowska-Mazur, K, & Lesch, H. 2004, ApJ 605, 
L33 

Hanasz, M., Kowal, G, Otmianowska-Mazur, K, & Lesch, H. 2006, AN 327, 
469 

Hawley, J.F, Gammie, C.F, Balbus, S.A. 1995, ApJ, 4 40, 442 

Hessen, V., Dettmar, R.-J., Krause, M., Beck R., 2008, arXiv:0801 .35421 
Jokipii, J.R.: 1999, in J. Franco and A. Carraminana (eds.) Interstellar 

Turbulence, Cambridge University Press, 70-78. 
Kowal, G., Hanasz, M., Otmianowska-Mazur, K, 2003, A&A, 404, 533 
Kowal, G, Otmianowska-Mazur, K, Hanasz, M., 2005, A&A, in press 
Lesch, H, Hanasz, M. 2003, A&A, 401, 809 
Maron, J., Blackman, E.G. 2002, ApJ, 566, L41 
Maron, J., Cowley, S., McWilliams, J. 2004, ApJ, 603, 569 
Otmianowska-Mazur, K, 2003, A&A, 408, 817 

Otmianowska-Mazur, K, Kowal, G., Hanasz, M., 2007, ApJ, 668, 1100 
Parker, E.N. 1992, ApJ, 401, 137 



